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Abstract 


The small polaron, a one-dimensional lattice model of interacting spinless fermions, 
with generic non-diagonal boundary terms is studied by the off-diagonal Bethe ansatz 
method. The presence of the Grassmann valued non-diagonal boundary fields gives 
rise to a typical U(l)-symmetry-broken fermionic model. The exact spectra of the 
Hamiltonian and the associated Bethe ansatz equations are derived by constructing an 
inhomogeneous T — Q relation. 
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1 Introduction 


In this paper we focus on constructing the Bethe ansatz solution of the small polaron with 
generic non-diagonal boundary terms, described by the Hamiltonian 


H 


^ ^— {cos{ri)nj+inj -f cos{r])nj+inj + c+cj+i d- c^+iCj] 
smr] 

+ - cot(-?/>_)[hi - ni] [k+un - H-Un] + csc{'iIj-) [a_ci /3-cf] 

-h csc{'ijj+) [a+CN + , (1.1) 


where cf and c,- are the creation and annihilation operators of spinless fermions at site j 
(which obey anticommutation relations {c^, Ck} = Sjk), respectively; the operators of particle 
numbers are nj = c^cj and Uj = 1 — nj; the parameters 'ijj±, a± and f3± are the boundary 
parameters related to boundary interactions; rj is the bulk coupling parameter. The boundary 
coupling K± is given by ^cscilj+cscrjsm^r] ± '0+) respectively. The model fll.ip is a typical 
spinless fermion model with boundary terms in condensed matter physics. It provides an 
effective description of the motion of an additional electron in a polar crystal DEI- In one 
spatial dimension, the model is integrable for both periodic and open boundary conditions 
by reconstructing it within the framework of quantum inverse scattering method (QISM) 

[alia El E]. 

In the past few decades, the integrability and the excitation spectrums problem have 
been studied extensively. For the small polaron model with periodic and purely diagonal 
boundary conditions, which can be mapped onto the XXZ quantum spin chain through the 
Jordan-Wigner transformation, the energy spectrum problem of the model was solved by the 
Algebra Bethe Ansatz method in [3181E]. A remarkable result was given by Yukiko Umeno [9] 
who constructed the fermionic R-operator and solved the spectrum problem via the Algebra 
Bethe Ansatz method. The generic integrable boundary conditions were obtained [10] by 
solving the graded reflection equation m- Subsequently, the Lax pair formulation of the 
generic integrable boundary conditions was presented in [12]. Since then, there have been 
numerous efforts to work out the exact solutions of the model. In 2013 the authors in [T3lfT^ 
figured out the Bethe ansatz solution of the model with non-diagonal boundary terms based 
on a deformation of the diagonal case and commented on the eigenstate of the model which 
envolves into the Fock vacuum when the off-diagonal boundary terms were ignored. The 
result is also closely related to that of algebra Bethe ansatz method. However, there still 
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exists a main obstacle for applying the conventional algebra Betlie Ansatz method to get the 
exact solution of the model with generic off-diagonal boundary conditions. The difficulty 
is mainly due to the fact the Hamiltonian fll.ip includes Grassmann valued non-diagonal 
boundary helds (or couplings) such as the terms associated with the parameters a± and 
which breaks the bulk [/(l)-symmetry of the model. The breaking of the [/(l)-symmetry 
leads to the obvious reference state (all-spin-up or all-spin-down state) is no longer the 
reference state in the usual algebraic Bethe ansatz [15] . 

Very recently, a systematic method for approaching the exact solutions of generic inte- 
grable models either with U{1) symmetry or not, i.e., the off-diagnonal Bethe ansatz (ODBA) 
method [16] was proposed in [13 UHl UHl |20] . With the ODBA method, some long-standing 
models [TH] |2T1 |22l [23l |2l] without U{1) symmetry were then solved. In this paper we study 
the small polaron model with the generic integrable boundary condition specihed by the 
X-matrices within Grassmann numbers via the ODBA method. 

The paper is organized as follows. In Section 2, we begin with a concise view of the 
integrability of the fermion model with the open boundary condition within the framework 
of the graded QISM. Some basic ingredients and algorithm of the transfer matrix are also 
introduced. In Section 3 we show that the Hamiltonian of the model can be rewritten 
in terms of the corresponding transfer matrix. In Section 4, after deriving the operator 
product identities of the transfer matrix at some special points of the spectrum parameter 
and its asymptotic behaviors, we express the eigenvalue of the transfer matrix in terms of 
an inhomogeneous T — Q relation and derive the associated Bethe ansatz equations. Finally, 
we summarize our results and give some discussions . 


2 Transfer matrix 


Let V be a two-dimensional Z 2 -graded vector space (or super space) [25] with an orthnormal 
basis = 1,2}. The grading of the basis vectors is [|1)] = 0, [|2)] = 1. The R-matrix of 
the small polaron model is given by [14] 


R{u) 


1 

sin 1 ] 


/ sin(M -|- r]) 


\ 


0 

0 

0 


0 0 0 \ 

sin u sin rj 0 

sin T] sin u 0 

0 0 — sin(u-l-p) y 


( 2 . 1 ) 
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acting on the tensor product V ®V oi two superspace. Here u is the spectral parameter and 
r] is the crossing parameter related to the bulk coupling fll.ip . The i?-matrix R{u) satishes 
the graded quantum Yang-Baxter equation (g-QYBE) [26] 


Ru{u - v)Ri3{u)R23{v) = R23{v)Ri3{u)Ri2{u - v), (2.2) 

and enjoys the properties: 

Initial condition : i?i2(0) = P 12 , (2.3) 

Unitarity relation : Ri 2 (u)R 2 i(-u) = ^(u), ^(u) = ^ ^ 2 . 4 ) 

sm r] sm rj 

P-symmetry : R2i{u) = Pi2Ri2{u)Pi2 = Ruiu), (2.5) 

T-symmetry : = R 2 i{u), (2.6) 

Crossing relation : R^fi—u — 2r])R2i{u) = ^{u + rj), (2.7) 

Antisymmetry : Rui—r]) = —2P^~\ (2.8) 

Periodicity : Ri 2 {u + vr) = —a-^i?i 2 ('u)ct^ = —a 2 Ri 2 {u)a 2 . (2.9) 


In the above equations, stj and istj are the partial super transposition and its inverse, Pij 
is the graded permutation operator and is a projector with rank one. 


p(-) 


/ 0 0 0 0 \ 

1 01-10 
2 0-110 
\ 0 0 0 0 / 


( 2 . 10 ) 


Here and below we adopt the standard notations: for any matrix A G End(Y), Aj is an 
embedding operator in the tensor space V ®V ® ■ ■ ■, which acts as A on the j-th space and 
as identity on the other factor spaces; Rij{u) is an embedding operator of i?-matrix in the 
tensor space, which acts as identity on the factor spaces except for the Tth and j-th ones. 
Since we discuss a fermionic lattice model, all embeddings are to be understood into a super 
tensor product structure. It is remarked that the super tensor product is graded according 
to the rule 


{A ® (2.11) 

where the parity [/)] is equal to zero (one) for bosonic (fermionic) indices. (For the details 
about the algorithm of super tensor product we refer the reader to [23 US]-) 
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We introduce two monodromy matrices To(m) and To('u), which can be considered as 2x2 
matrices on the auxiliary space with elements being operators acting on , 

To{u) = Ron{u — 0]s[)Ron-i{u — O^^i) ■ ■ ■ -Roi('^ — ^i), (2-12) 

To{u) = Roi{u + 9i)RQ2iu + O 2 ) ■ ■ ■ Ron{u + 0^). (2-13) 

Here {Oj\j = 1, 2, • • •, A^} are arbitrary free complex parameters which are usually called the 
inhomogeneous parameters. 

The framework of QISM for integrable systems with open boundary conditions in a way 
that makes it applicable to super spin chains. Following EE], for a given R-matrix, we 
introduce a pair of K-matrices K~{u) and K^{u). The former satishes the graded reflection 
equation 

Ri 2 {u - v)K{{ u)R 2 i{u + v)K 2 (v) 

= K2{v)Ri2{u + v)K^{u)R2i{u-v), (2.14) 

and the latter satishes the dual graded rehection equation 

Ru{v - u)Kt{u)R 2 i{-u - 

= KY{v)Ri2i-u - vY^^^'^^^KY{u)R2i{v - u), (2.15) 

whereas the new matrices R and R are related to the R-matrix via 

4 ,( 2 . 16 ) 

^^^(^)isii,.i2 _ ([{Rr2^(w)K*']"^)***h (2.17) 

For open super spin chains, rather than the standard monodromy matrix To('u) fl2.12p . 
we need to consider the double-row monodromy matrix To(-u) 

To{u) =To{u)KQ{u)fo{u). (2.18) 

Then the double-row transfer matrix t{u) of the system is given by 

t{u) = sfro{iF(^(w)To(M)}, (2.19) 

where str{-} denotes the super trace of a matrix, which is dehned by 

str{A} ^ 5^(-1)WH). (2.20) 
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The graded QYBE fl2.2l) and REs fl2.14p and fl2.15p lead to the fact that the transfer matrices 
give rise to a family of commnting operators |6] with different spectral parameters: 

[t{u),t{v)] = 0. (2.21) 

Then t{u) serves as the generating fnnction of the conserved qnantities, which ensnres the 
integrability of the system. 


3 Small polaron with open boundaries 


In this paper, we consider the K-matrices K (u) and K~^{u) which satisfy the graded REs 
mu and possess the following generic expressions (see also uniEHipig) 


K-(u) = u- 


sin (« + ?/)_) Q;_sin(2'u) 

/5_sin(2n) — sin(n — 


K~^{u) = UJ+ 


sm{u + r] + a_|_ sin(2[n + 77]) 
sin(2[M + rj]) sin(M + r] — 

— 1 


(3.1) 


(3.2) 


with normalizations u± dehned by ^ and 0J+{ri) = 2cos(r;/sin(p4.) • Here a±, 

/9± are all Grassmann nnmbers which are related to bonndary fields. The parameters 'ip± are 
arbitrary commnting even Grassmann nnmbers bnt the invertibility reqnires them to have 
a non-vanishing complex part, the remaining non-diagonal bonndary parameters a± and /3± 
are anticommnting odd Grassmann nnmbers, namely. 


IV’+i '/'-]= 0 = {a±, = {a±, &} = {/?±, /3±}. 


(3.3) 


In addition, the odd Grassmann nnmbers are snbject to the condition a±f3± = 0 dne to the 
graded REs fl2.14l) and fl2.15l) . 

Based on the graded QISM, the Hamiltonian fll.ll) of the small polaron model with generic 
off-diagonal bonndary terms can be rewritten in terms of the transfer matrix fl2.19p as: 


H = 


1 dt{u) 1 

^^\u=o,{e,=o} + 


N-l 


stro{KQ (0)} -I- ^ stro{R'(('(0)PAroRoAi(0)} 

j=i 


+K^ (0 )-h-tan? 7 . 


(3.4) 


The pnrpose of this paper is to constrnct the spectra of the Hamiltonian and derive the 
corresponding Bethe ansatz eqnations. 
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4 Eigenvalues and the Bethe ansatz equations 

4.1 Functional relations 

Following the similar method developed in |20], we derive that the prodncts of the transfer 
matrix fl2.19p of the snper spin chain with the generic open bonndaries described by the 
K-matrices in fl3.ip and fl3.2l) . at the points 6j and 6j — rj, satishes the relations 

W - ^ = 1,...^. (4.1) 

For generic {^j}, the qnantnm determinant operator Ag{6j) is proportional to the identity 
operator, namely, 


Aq('u) = 6{u) X id, 


(4.2) 


where the fnnction 6{u) is given by 




5{u) = —2— sin(M + '0+) sin(n — sin(M + sin(n — 'ijjJ) sin(2n + 2rj) sin(2n — 2r]) 
sin rj 


X 


N 

n 

1=1 


sin('u — 9i — rf) sin(n — 6i + rf) sin(M + 6i — rf) sin('u + 6i + rf) 


sin^ r] 


sin^ rj 


(4.3) 


Fnrthermore, we have checked that the transfer matrix t{u) of the small polaron model with 
the generic bonndary conditions enjoys the crossing property 

t{—u — rj) = t{u). (4.4) 

The qnasi-periodicity of the R-matrix (12.91) and K-matrices 

i?i2(n -|- tt) = —alRi 2 {u)(^i = —(^ 2 ^i 2 iu)a 2 ', K^(u -|- vr) = —a^K^(u)a^, (4.5) 

and the special points valnes at n = 0, | of the K-matrix give rise to several properties of 
the associated transfer matrix, namely. 


t{u + T{) = t{u), 


N 


*(o) = n 


sm{r] - 6i) sin(r^ 9i) 


1=1 


sin^ rj 

N 


X id, 


,x/l, , , ,, -p-rsin(f-0; + 77)sin(f+ 0; + //) 

t{—) = cot yi-cot- - - — - - -x id. 


i=i 


sin^'/; 


(4.6) 

(4.7) 

(4.8) 


.lim tin) = X U1 (4.9) 

iu=i-±oo [21)"^^'^^ sin Tj 
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Here the operator is given by 


N 

i=i 


a. 


3 ’ 


(f/*)" = id, 


(4.10) 


which commutes with the transfer matrix. The relation fl4.10p allows us to decompose the 
whole Hilbert space "H into two subspaces, i.e., "H = ®'H~ according to the action of the 

operator U^\ The commutativity of the transfer matrix and the operator 

U^, i.e., [t{u), U^] = 0, implies that each of the subspace is invariant under t{u). Hence the 
whole set of eigenvalues of the transfer matrix can be decompose into two series, denoted 
by A±(ti) respectively. The eigenstates corresponding to A+(ri) (resp. A_{u) ) belong to the 
subspace (resp. "H”). The operator product identities fl4.ip of the transfer matrix and 
the commutativity of the transfer matrix with different spectrum u enable us to derive the 
following relations of the associated eigenvalues A±{u) respectively. 


A±{0j)A±{ej - r]) 


6{6j) smrjsmr] . ^ ^ 

sm{29j + rj) sin{29j — r])' ^ 


(4.11) 


with the function S{u) given in fl4.3p . 

The properties of the transfer matrix t{u) given by fl4.4l) - fl4.9p . imply that the correspond¬ 
ing eigenvalue functions A±(m) satisfy the relations: 


A±{-u - rj) 


A±(0) 



lim A±(m) 


A±(m), A±(M-h7r) = A±(m), 


N 

n 

1=1 


sm{r] — 9i) sin(p -|- 9i) 


sin^ r] 


N 


cot l/j- cot -0+ 


sin(| - 9i + r]) sin(| + 9i + rj) 


1=1 


sin^ rj 


± U3j^U3- 


(2z; 


2A1+2 


(a+/3_ - P+a_) . 


2N 
Sin T] 


(4.12) 

(4.13) 

(4.14) 

(4.15) 


Obviously, A±(w) are a degree 2A^-|-4 trigonometric polynomial of u, along with the crossing 
symmetry and the periodicity mentioned in fl4.12p . these factors lead to that only A^ -|- 3 
unknown coefficients need to be determined by -|- 3 special points values of the associated 
function A±(m). Therefore, the two functions A±(m) can be completely determined by the 
above functional relations fl4.lip - fl4.15p . 





















4.2 Eigenvalues of the transfer matrix 


Following the method in [T^ HU HHl |20] and with the helps of the fnnctional relations fld.lip - 
fl4.15p . we can express the eigenvalne A±(ti) of the transfer matrix of the small polaron model 
with the bonndary terms specihed by the generic non-diagonal K-matrices given by fl3.ip and 
fl3.2p in terms of an inhomogeneons T — Q relation [18] respectively, 


A ^ ^ ^ ^Q^^Ku-t]) ^Q^^\u + r]) 

A±(n) = a{u) — , ' -|-a(n)- 


± 


csin( 2 M) sin( 2 M -|- 2r]) 




A{u)A{—u — r]), 


(4.16) 


where the Q-fnnctions are parameterized by | j = 1, • • •, A^} respectively 

sin('U + Tj) 


Q“(«) = If 

1=1 


sm m 


= - 7]). 


sm Tj sm T] 

The other fnnctions a(n), d(n), A{u) and the constant c are given by 


(4.17) 


N 


A{u) = n 

l=l 


sin(M — 9i + 7]) sin(n -|- 6 *; -|- 77 ) 


sin^ 1] 


7 \ ■ ( I \ ( I 4sin(2n + 2r7) ^ 

a{u) = uj+u}-sm[u — 7p+)sm[u — ip-) - - -— Apn), 

sin( 2 'u -|- 7]) 

7/4 • / , 4 . / ,4 sin( 2 M) - 

d[u) = uj+U-Sin[u + T] + Tp^) sm[u + T] + — ^—Ay—u — rj) 


sin( 2 n -|- 77 ) 


= a(-M-77), 

c = a;+a;_(Q(+/5_ —/9+a_). 


(4.18) 

(4.19) 

(4.20) 

(4.21) 


Since that A±(m) both are polynomials, the residnes of A±(n) at the apparent poles u = 
and u = —~ V, j = '' T ^ mnst vanish, which gives rise to the associated BAEs 

«(Af ^ -V) + ^ + V) 

±csin2/7^."^^ sin(2/7^.'^^ 2T])A{fif"^)A{-ij.f'^ - 77 ) = 0, j = 1,... ,N. (4.22) 

It is easy to check that the T — Q relation fl4.16p satisfies the relations fl4.11l) - fl4.15l) re¬ 
spectively nnder the condition of N parameters {fij \ j = 1,---,A^} satisfying the BAEs 

(1122D. 


9 

























We remark that the roots = 1, • • •, N} to the BAEs (I4.22p are Grassmann number 


valued, which implies that the corresponding Q-functions in fl4.17p can be expressed as 

Q^^\u) = Qit\u) +gQ^^\u), g = a+j3_ - /3+a_, andgf^ = 0, (4.23) 


where 


of’ 


sin(u - sin(M + Af 

aW) = 11- 

i=i 


smr; 


smr; 


N-l 


QfHn) = A^'II 


sin('u — sin(n + + g) 


i=i 


sm g 


smg 


(4.24) 


The 2N parameters = 0,1; j = are c-number valued. Substituting the 

relations fl4.23p - fl4.24p into the BAEs fl4.22p . one may get the associated 2N BAEs. The re¬ 
sulting BAEs completely determine the 2N c-number valued parameters = 0,1; j = 

1 ■ ■ ■ ,N}, which resembles those in [T3l [14] . 

In the homogeneous limit 6j —)■ 0, the above BAEs become 

+g)\ - -0+) - ?/)_) sin(2pj"^^ -F 2g) 


T- 


_ Q(/if ^ + g) 

^ + g ++ g + ^jJ_)sm{2|J.f^^) - v) 

(a+/5 - /?+«-) sin(2;Uj'^^ -h g) + 2g) sin^^(/ij'‘'^ -h g) 


sm /ij 


sin(/x^"^^ + g + '0+) + g + ip-) sin^^ gQ{gf^^ - g) 


j = l,...,iV. (4.25) 


Then two series eigenvalues of the Hamiltonian fll.ip can be expressed in terms of the Bethe 
roots as follows 


E+ = 


1 cIA-i-(m) 


u=o,{6»j=o} +-tan77 


2 du 
= -^cotip+- - cot ip- 


-h N cot 7/ - tan 7] 

sin 2g 2 


N 


sin^ g 


sin gj ^ sin(;uj + g) 
where the parameters = 1,..., N} satisfy the associated BAEs fl4.25p . 


(4.26) 


5 Conclusions 


The small polaron model with off-diagonal boundary described by the K-matrices fl3.1l) and 
h^ . which can be regarded as a graded version of the general open XXZ spin chain, has 
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been studied by the off-diagonal Bethe ansatz method proposed in [TH [T^ [HI [20] . Based 
on some properties of the R-matrix and K-matrices, we obtain the operator identities f|4.ip 
of the transfer matrix and then construct the corresponding inhomogeneous T — Q relation 
for its eigenvalues fl4.16p and the corresponding BAEs fl4.22p . Moreover, the exact spectra 
of the Hamiltonian is given in fl4.26p . When the nondiagonal boundary parameters satisfy 
the constraint a± = (3± = 0, the resulting T — Q relation is reduced to the conventional one 
which is the solution of the model with diagonal boundaries. 

A possible extension of the present work is to consider the mulit-component Bose-Fermi 
mixtures with off-diagonal boundary conditions with the help of the fusion method |29j . 
Meanwhile, according to the spin-s XXZ Heisenberg chain with generic non-diagonal bound¬ 
aries solved in [30], the construction and the solution of graded higher spin chain may be 
obtained by similar method. 
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